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ABSTRACT 

So far essentially all black hole masses in X-ray binaries have been obtained by 
observing the companion star's velocity and light curves as functions of the orbital 
phase. However a major uncertainty is the estimate of the orbital inclination angle of 
an X-ray binary. Here we suggest to measure the black hole mass in an X-ray binary 
by measuring directly the black hole's orbital motion, thus obtaining the companion 
to black hole mass ratio. In this method we assume that accretion disk wind moves 
with the black hole and thus the black hole's orbital motion can be obtained from 
the Doppler velocity of the absorption lines produced in the accretion disk wind. We 
validate this method by analyzing the Chandra/HETG observations of GRO J1655-40, 
in which the black hole orbital motion (Kbh = 90.8 ± 11.3 km s _1 ) inferred from the 
Doppler velocity of disk-wind absorption lines is consistent with the prediction from 
its previously measured system parameters. We thus estimate its black hole mass 
(M BH = 5.41±oj7 M©) and then its system inclination (i = 72.0t 7 7 8 5 °), where M BH 
does not depend on i. Additional observations of this source covering more orbital 
phases can improve estimates on its system parameters substantially. We then apply 
the method to the black hole X-ray binary LMC X-3 observed with HST/COS near 
orbital phase 0.75. We find that the disk-wind absorption lines of C iv doublet were 
shifted to ~ 50 km s _1 , which yields a companion-to-black-hole mass ratio of 0.6 
for an assumed disk wind velocity of —400 km s _1 . Additional observations covering 
other orbital phases (0.25 in particular) are crucial to ease this assumption and then 
to directly constrain the mass ratio. This method in principle can also be applied 
to any accreting compact objects with detectable accretion disk wind absorption line 
features. 



1 INTRODUCTION 

Black holes (BHs) are believed to exist in many X-ray bi- 
naries (XRBs) and active galactic nuclei (AGN). Measure- 
ment of the motion of a BH with respect to its surrounding 
in such a system can in principle set strong constraints on 
the mass of the BH. However in an AGN this is normally 
not possible because the time scale of BH's significant mo- 
tion is too long and/or the BH is barely moving at all. In 
an XRB, the BH should move with respect to the system's 
center of mass (CM), making it possible to detect directly 
the BH's motion with respect to the CM. However in a low- 
mass XRB (LMXB) the companion is normally much less 
massive than the BH, so that the BH barely moves or moves 
very slowly with respect to the CM. In some LMXBs, the 
companion's mass is comparable or only sever al times less 
mass i ve than the BH, e.g., in GRO J1655-40 (|Zhang et al.l 
1 19941 ; lOrosz fc Bailvnlll997lh the BHs' motion may be sig- 



nificant enough for direct detection. The most favorable sys- 
tems for detecting BH's motion should be high-mass XRBs 
(HMXBs), in which the BHs move rapidly with respect to 
their CMs. 

Orbital motion of double-peaked disk emission 
line s were observed for n eutron star XRB Sco X- 
1 ilSteeghs fc Casaxesl l2002r), the BH X RB A0620-00 
l|Haswell fc Shafterl Il990l: lOrosz et all Il994l ). and the BH 
GRS 1124-68 |Orosz et al.lll994h .~Unfortunatelv a signif- 
icant phase offset of velocity modulation was found from 
that expected based on the observed orbital motion of the 
companion, though the veloci ty semi-amplitude is consistent 
with the expected mass ratio (|Orosz et al.lll994h . Soria et al. 
(1998) observed the orbital motion of the double-peaked disk 
emission line He n A4686 from GRO J1655-40, and found 
its velocity modulation phase and semi-amplitude in agree- 
ment with the kinematic and dynamical parameters of the 
system. Therefore a more robust mass lower limit is placed 
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based on the observed motion of the primary and thus ruling 
out any possibility for a neutron star as the primary in the 
system (Soria et al. 1998). However one major problem in 
accurately measuring the orbital motion of the primary from 
the observed double-peaked emission lines is how to deter- 
mine reliably the line center, because the lines are typically 
asymmetric and also variable. 

In an XRB, both the accretion disk and its wind move 
with the BH, and thus provides us with another opportunity 
to measure the BH's motion via Doppler shift of absorption 
features of the accretion disk wind. Accretion disk winds are 
ubiquitous in XRBs and normally detected through ionized 
absor ption lines, typically with around ~ 10 3 km s -1 o r less 
(e.g.. lUeda etld1l2004l ; [Miller et alJliooj , 120061 . liooj ). but 



can reach to about 



10 km s in some extreme cases, 



e.g. i n the newly dis covered BH transient IGR J17091- 
3624 l|King et alj 201 ll ) . In particular the high quality Chan- 
dra/HETG observations of the BHXB GRO J 1655-40 have 
found many highly ionized narrow absorption lines, which 
are interpreted as evidence of magnetic field-driven ac- 
cretion disk wind dMiller et all 120061 ; iKallman et"aH 120091 ; 
iLuketic et al.|[2010l ); although the absorption lines are also 
interpreted as from the absorption by X-ray heated thermal 
wind (Netzer 2006). Regardless the origin of the accretion 
disk wind, the orbital motion of its many absorption lines 
may be measured reliably, because typically many narrow 
absorption lines are present with high signal to noise ratios 
and appear to be rather stable when observed. 

In this work, we model the Doppler motion of wind ab- 
sorption lines from the LMXB GRO J1655-40 and HMXB 
LMC X-3 to constrain directly the companion to primary 
mass ratio, in order to measure their BH masses and or- 
bital inclination angles. We first describe the methodology 
and test its feasibility by applying it to Chandra/HETG 
observations of GRO J1655-40, revealing for the first time 
the velocity modulation of wind absorption lines in an XRB 
and thus providing a new measurement of its BH motion 
and orbital inclination angle. We then apply this method 
to the HST/COS observations of LMC X-3 attempting to 
constrain the companion to primary mass ratio. Finally we 
discuss further observations needed to achieve the required 
accuracy of BH mass estimate for LMC X-3, as well as po- 
tential problems and uncertainties of applying this method. 



2 METHODOLOGY 

So far, all BH masses in XRBs have been estimated using 
the Kepler's 3rd law of stellar motion, expressed in the so 
called the mass function, 



f(M) EE P olh K c /27vG 



M BH sin 3 i/(l + g) 2 , (1) 



where P or b is the orbital period, Kc is the semi-amplitude of 
the velocity curve of the companion star, Mbh is BH mass, 
i is the the orbital inclination angle, and q = Mc/Mbh 
is the mass ratio. Since the only direct observables are P OT b 
and Kc, both Mc and i have to be determined indirectly in 
order to obtain the BH mass estimate reliably. The compan- 
ion's mass Mc can be determined relatively reliably by the 
observed spectral type of the companion star. For LMXBs, 
i can be estimated by modeling the optical or infrared light 



curve modulation, though model dependence and other un- 
certainties (such as accretion disk contamination) cannot be 
circumvented completely. For HXMBs, i is normally not de- 
termined very well; in many cases observations or lack of 
eclipse of the accretion disk emission by the companion is 
used to put some constraints on the possible ranges of i. For 
de tails of BH mass estimate s using this method, please refer 
to lRemillard fe McClintockl i|2006l i. 

On the other hand, the mass ratio q can be determined 
directly according to the law of momentum conservation, 



Mc/Mbh = Kbb/K c , 



(2) 



if the semi-amplitude of the velocity curve of the BH Kbh 
can be observed directly. Since a BH is not directly observ- 
able, we can only hope to observe any emission or absorption 
line feature co-moving with it. The accretion disk certainly 
moves with the accreting BH. However any line feature of 
the inner accretion disk suffers from the broadening of disk's 
orbital motion and distortions by relativistic effects around 
the BH, thus making it practically impossible, or difficult to 
the least, for detecting the binary orbital motion of the BH. 
Emission line features from the outer disk region are nor- 
mally detected with double-peaks, which can be modeled to 
obtain the semi-amplitude of the velocity of the compact 
object, as discussed above. However a major uncertainty is 
to determine the mean separation between the emission re- 
gions of the blue-shifted and red-shifted components. 

Fortunately, as we have discussed above, absorption line 
features of accretion disk winds in BHXBs have been rou- 
tinely detected with high significance. This suggests that 
the accretion disk wind in an XRB moves with the disk that 
produces the wind, since otherwise the wind from the disk 
would not intercept the continuum emissions produced from 
the same disk. In this case, the Doppler motion of wind 
absorption line features can be considered as that of the 
BH, unless the wind interacts strongly with the surrounding 
interstellar medium or the wind intrinsic velocity has sys- 
tematic orbital dependence. It has been found recently that 
large scale (pc) cavities exist around microquasars (XRBs 
producing relativistic jets), and perh aps are ubiquitous i n 
all BHXBs producing strong winds (|Hao fc Zhand |2009h . 
This suggests that at least in the vicinity of the BH the 
wind cannot interact directly with interstellar medium. The 
interactions of disk winds with interstellar medium at pc 
scale would not complicate the observations, because it takes 
much longer than an orbital period before the winds could 
arrive at the boundaries of these cavities, and thus the wind 
should have lost any possible memories of the orbital motion 
of the compact object. For disk-fed accretion, the wind ve- 
locity is not expected to have any orbital phase modulation; 
for stellar-wind accretion the focused wind may complicate 
the situation considerably though. We therefore suggest to 
apply this method to only those systems in which the accre- 
tion is disk-fed, i.e., no significant wind interaction happens 
between the stellar and accretion disk winds. Nevertheless 
only observations can tell us in which kinds of systems this 
method can be used reliably. 
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Figure 1. Distributions of the line velocities in four equally 
spaced intervals and their average values, during the Chan- 
dra/HETG observations lasting for more than 60 ksec; note that 
here the system velocity of the BH XRB has already been sub- 
tracted. Panels from (a) to (d) corresponding to the phases shown 
in Fig.[2]from left to right. 
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Figure 2. Velocity curve of the 39 observed absorption lines, 
after rejecting six lines as suspected outliers and subtracting the 
line of sight intrinsic velocity of each line at each orbital phase. 
The upper panel marks the velocity of each line with its 1-rr error 
bar slightly shifted horizontally for visual clarity; the inset shows 
all velocities, including the three data points out of the range in 
the main panel. The bottom panel shows the weighted average 
velocity of all lines in the upper panel at each phase; the solid 
curve is the fitted velocity curve with its orbital period and phase 
fixed at the values observed previously. 



3 APPLICATION TO THE LMXB 
GRO J 1655-40: A TEST STUDY 

GRO J1655-40, discovered as an X-ray transient by 
Zhang et al. (1994), is a well-known and best studied 
BHXB and the se cond source with superluminal relativis- 
tic jets detected dTingav et al.l 1 19951 ; lHarmon et al. I Il995l ; 
iHiellming fc Rupenl 1 19951 ). Its system parameters remain 
so far the best measured among all known BHXBs, with 
P orb = 2.62191 ± 0.00020 d, 1/q = M BH /M C = 2.6 ± 0.3, 
70.2±1.9°, and M B h = 6.3±0.5 Mp (al l 95% confidence) 



i|Orosz fc Bailvnlll997l ; iGreene et al.ll200ll ). Its precise BH 
mass and inclination measurements allow its BH spin pa- 
rameter determined from its X-ray continuum fitting, first 
proposed by Zhang et al. (1997) and then refine d by incor- 
porat i ng detailed modeling of various effects (e.g..lYao et al 
20051; iLi et "ail 120051; iDavis et alj 120051; IShafee et al.l I200T 



Steiner et al.l 120091 ; (McClintock et al.l l201l[ ). However pos- 



sible sources of systematic errors are in the modeling of the 
ellipsoidal light modulation of the companion star and con- 
tamination of the disk's continuum emission in the optical 
to infrared bands. Soria et al. (1998) observed the velocity 
curve of the double-peaked disk emission line He 11 A4686 
from GRO J1655-40, and derived its projected radial veloc- 
ity semi-amplitude determined for the primary as Kbh = 
76.2 ± 7.5 km s~\ yielding M BH = 6.62 ± 0.74 M© and 
i = 66.6 ± 7.7° (all 68.3% confidence), fully consistent with 
that determined from the companion's velocity curve and 
ellipsoidal light variation. 

Miller et al. (2006) detected many narrow and ionized 
absorption lines with High Energy Transmission Grating 
(HETG) Chandra observations of GRO J1655-40. The exact 
position of the point source cannot be determined accurately 
in the zeroth order image, which is severely piled-up due to 
its high flux. An offset in the source position will cause a 
systematic shift of all lines in the wavelength space. How- 
ever, this shift is mostly canceled out in the combined total 
spectrum if each line is detected in the spectra of both the 
+ 1 and —1 orders with similar count spectra. However the 
combined lines will be broadened by the unknown system- 
atic offset, which in turn will cause larger uncertainties in 
determining the center of each line. There may also exist 
other sources of systematic errors which can cause offset be- 
tween the centers of the same line in the spectra of the +1 
and —1 orders. 

Here we determine the position of the zeroth order im- 
age by minimizing the relative offsets between the absorp- 
tion lines in the +1 and —1 orders, and estimate the sys- 
tematic error in determining the center wavelength of each 
absorption line. We first fit each absorption line feature with 
a Gaussian profile in the +1 and —1 orders of the High En- 
ergy Grating (HEG) and Medium Energy Grating (MEG) 
separately, i.e., we get a pair of profiles for each absorption 
line, integrated over the whole observation in order to get 
the highest signal to noise ratio. The local continuum around 
each line is assumed to have a power-law shape, but our re- 
sults are insensitive to the continuum shape. For HEG and 
MEG we obtain 43 and 9 pairs of absorption lines with at 
least 3-cr detection, resulting in 43 and 9 differences of the 
central wavelengths of these pairs, respectively. We then use 
the Bayesian method to find the offset d and its systematic 
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error ct s by maximizing the following function, 



p(d,a B ) = 



1 



n 



(3) 



where Xk and o Xk are the wavelength difference and its sta- 
tistical error of the fc-th pair (k = 1 to n, and n = 43 for 
HEG or 9 for MEG), respectively. Finally for HEG we get 
d = (4.8 ± 0.1) x 1CT 3 A and a s = (3.2 ± 1.3) x 1CT 4 A. 

We thus take the error for d as ad — cr B — 3.2 x 1CP 4 A 
for HEG. Similarly for MEG we get d = (9.7 ± 0.8) x 10~ 3 
A for MEG, where <7<j = 9 x 10~ 4 A is dominated by the 
systematic error. Taking these additional errors for MEG 
and HEG line wavelengths and requiring all pairs in HEG 
with d = 4.8 x 10" 3 A and in MEG d = 9.7 x 10" 3 A, 
we have a total \ 2 ~ 55 for 50 degrees of freedom, thus 
validating our method of determining the systematic errors. 
Since the wavelength bin (an image pixel) is 2.5 x 10~ 3 A or 
5.0 x 10~ 3 A for HEG or MEG, we thus shift the spectra of 
the —1 and the +1 orders by one bin, towards the decreasing 
(— 1 order) and increasing (+1 order) wavelength directions, 
respectively. This is identical to shifting the location of the 
point source by one pixel in the zeroth order image. Finally 
we combine the two shifted spectra as one single spectrum. 

Because the observations lasted for about a quarter of 
the orbital phase of GRO J1655-40 and many absorption 
lines have high signal to noise ratios, here we divide the ob- 
servations into four equal orbital internals in order to detect 
the orbital motion of the absorption lines. For each com- 
bined spectrum we again fit each absorption feature with 
a Gaussian profile to determine its central wavelength, line 
width and intensity. Table 1 lists the 45 lines in all four 
intervals with at least 3-<x detection that are included for 
further analysis; here only statistical errors are shown. 

Fig.Q]shows the velocity distributions of the four groups 
of 45 lines and their average velocities; note that here the 
system velocity of the BH XRB has already been sub- 
tracted. Adopting the orbital period and phase ephemeras 
from Greene et al. (2001) and assuming that the BH's mo- 
tion is exactly anti-phased with the companion's motion 
(Soria et al. 1998), we fit velocities of these lines simul- 
taneously to a sinusoidal function. In the fit, the intrinsic 
velocity of each line is a free parameter, but all lines in 
the same orbital phase follow the same orbital modulation. 
This way we have 46 free parameters (45 intrinsic velocities 
for all these lines plus the line of sight velocity amplitude 
of the black hole) with 180 data points. The fit results in 
Kbh = 93.8 ± 11.1 km s _1 with x 2 = 179 for 134 degrees 
of freedom; the systematic errors determined above are in- 
cluded for all lines. 

Considering that some lines may be outliers, we reject 
the four groups of six lines with central velocities more than 
300 km s _1 from the median values of each distribution, 
as listed in Table 1. The remaining four groups of 39 lines 
have velocity dispersion (1-cr) of about 60 km s _1 , i.e., the 
rejected lines are more than 5-<r away from the median val- 
ues. The fit to the remaining four groups of 39 lines yields 
Kbh = 90.8 ± 11.3 km s _1 with \ 2 = 131 for 116 degrees of 
freedom, a marginal improvement over the full dataset fit. 
This means that the deviations of these possible outliers are 
not significant statistically and thus do not deserve further 
studies at this stage. Fig. [2] shows the fitting results; please 
note that the velocities are obtained by subtracting the fitted 
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Figure 3. Spectrum of the available HST/COS observations 
around C IV doublet. Red curves indicate the fit to the emis- 
sion from the heated surface of the companion star, green curves 
are the fit to the Galactic ISM absorptions at km s —1 , and blue 
curves indicate the fit to the accretion disk wind absorptions at 
50 km s~ x (local velocity S3 —400 km s _1 ). The vertical bars 
mark the central positions of these components. 



line of sight intrinsic velocity of each line from the obtained 
central velocity of each line at this orbital phase. It is worth 
noting that the phase zero in this system was defined as that 
when the companion is receding from the observer at the 
maximum velocity (Orosz & Bailyn 1997). Nevertheless both 
results are statistically consistent with Kbh = 76.2 ±7.5 km 
s" 1 obtained by Soria et al. (1998). If we fit the four groups 
of 39 lines to a straight line, we get \ 2 = 813 for 154 de- 
grees of freedom. The linear model is thus rejected with high 
significance, compared to the sinusoidal model. 

With the fitted Kbh = 90.8 ± 11.3 km s" 1 and tak- 
ing the other system parameters (except the inclination) of 
GRO J1655-40 from Greene et al. (2001), we first obtain 
its BH mass M BH = 5Al±°- 9 £ M Q from Eq. 2 and then 
its system inclination i — 72.OI75 from Eq. [1] respec- 
tively. All these (l-o") errors are obtained by Monte-Carlo 
samplings, because of the asymmetry and coupling of some 
errors. These parameters, albeit with large uncertainties due 
to the very incomplete orbital coverage of observations, are 
consistent with all previous measurements. Therefore the 
existing data of GRO J 1655-40 validate our proposal that 
absorption lines produced in the accretion disk wind can be 
used to measure directly the orbital motion of the BHs in 
BHXBs. 



4 APPLICATION TO THE HMXB LMC X-3 

LMC X-3 is another excellent object for applying this 
method. It is a bright XRB system in the Large Magellanic 
Cloud composed of a B3 V star and a central BH, and reced- 
ing away from us at a systemic velocity Vs = +310 km s _1 
l|Cowlev et al. I Il983l ; hereafter C83). It is one of the few 
BH systems that are persistently luminous in both X-ray 
and far-ultraviolet (FUV) wavelength bands. The BH is be- 
lieved to be undergoing accretion from its B-star compan- 
ion via Roche lobe overflow with an orbital period of 1.7 
days (C83); we therefore do not expect any significant con- 
tamination of the stellar wind to any absorption features 
of its accretion disk wind. Spectroscopic observations of 
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Figure 4. Si II and Si IV absorption features in the spectrum 
of the previous HST/COS G130M observations, normalized to 
the continuum flux. The apparent higher velocity components 
(marked with pink curves) may be produced in the outer part of 
the disk wind at lower local velocities. 

the B star indicate a large radial velocity semi-amplitude, 
Kc = 235 ± 11 km s" 1 (C83). Taken the mass of the B3 V 
star as about 6 Mq, the BH mass in LMC X-3 is thus es- 
timated to be 5 — 10 Mp,, assuming an inclination an gle of 
50° - 70° (C83; iKuiper et all [l98s[ ISoria et al]|200ll ). Be- 
cause of the considerable uncertainty in its inclination (and 
thus BH mass), its BH spin has not been reliably deter- 
mined yet with the X-ray continuum m odeling method (e.g., 
IZhang et aflll997l : lsteiner et alJl20ld ). Since the BH should 
be moving at a comparable velocity to its companion, it is 
possible to directly measure its orbital motion if the absorp- 
tion lines produced in its accretion disk wind are detected. 

Although LMC X-3 has been intensively studied in the 
X-ray band, the previously operated and currently operat- 
ing X-ray spectroscopic instruments lack the combination 
of the required sensitivity and resolution to measure the 
expected Doppler motion of any accretion wind absorption 
features. High resolution UV spectroscopic observations of 
this source are rather sparse; so far only two FUSE ob- 
servations and four observations with the Cosmic Origins 
Spectrograph (COS) onboard the Hubble Space Telescope 
(HST) (HST/COS her e after) to LMC X-3 have been made 
jHutchings et al.ll2003l ; IWang et al.ll2005l ; ISong et allkoid . 
here after S10). Again FUSE does not have the combination 
of the required sensitivity and the observations did not pro- 
duce spectra with high enough signal to noise ratio for the 
purpose of this investigation. The HST/COS does have the 
required performance to do so, although the available obser- 
vations only covered a small portion of the orbital period of 
LMC X-3. In this work, we analyze the available HST/COS 
UV spectroscopic observations of LMC X-3, aiming at con- 
straining the systemic parameters and further demonstrat- 
ing the feasibility of our proposed new method of measuring 
the BH's mass in an XRB. 

From the FUSE and HST/COS observations, S10 de- 
tected variable O VI and N v emission lines and found that 
variability of their intensities are inconsistent with expecta- 
tion of a stellar wind origin. So they attributed them to the 
heated stellar surface of the companion star. In fact both 
absorption and emission features are detected in N v (with 
G130M) and C iv (with G160M) doublets (Fig. 2 in S10) in 
the HST/COS observations (Program 11642; please refer to 



S10 for detailed description of the HST/COS observations 
and data analysis). However, since S10 focused on the ori- 
gin of the emission variability and the connection between 
the emission features and the systemic parameters, these 
absorption features have not been explored in detail. 

Here we analyze the C iv doublet complexes (we do not 
conduct the same analysis to the N v features because the 
absorptions are much less significant than those of C iv). 
In Fig. [31 we decompose the observed C iv doublet into 
three components for each line: one emission component and 
two absorption components. The emission component is as- 
sumed to come from the heated surface of the companion, 
the same as the O VI and N v emission lines (S10); its width 
and Doppler shift are fixed to those inferred from the O VI 
and N v emission lines in the decomposition. One absorp- 
tion component (at zero velocity) is assumed to come from 
the Galactic ISM absorption. The other significant absorp- 
tion component has a redshift of about 50 km s _1 , and may 
have two possible origins. One is the absorption produced in 
the stellar wind of the B3 V companion star. However, the 
orbital phase of the binary system was about 0.75 during 
the HST/COS observations, i.e., the straight line connect- 
ing the companion star and the BH is perpendicular to our 
LOS (please be noted that, in contrast to GRO J1655-40, 
the orbital phase zero in LMC X-3 was defined as that when 
the BH is at its superior conjunction, i.e., the companion is 
just between the observer and the BH; Cowley 1983; S10). 
In order for any stellar wind to intercept significantly the 
emission from the inner disk region, a significant amount 
of stellar wind has to stream to the inner disk region, i.e., 
the system must be wind-fed, against the common believe 
that LMC X-3 is actually disk-fed. We thus rule out this 
possibility. The other and the only viable scenario is the ab- 
sorption by the accretion disk wind, making it possible to 
directly measure the orbital motion of the BH. It is worth 
noting that the measured C iv column density in the ISM 
[log A r isM(cm -2 ) = 13.7 ± 0.1] agrees remarkably well with 
predicted value [log JVciv(cm -2 ) = 13.8] in a joint analysis 
of the X-ray and FU V spectroscopic observations of LMC X- 
3 (|Yao et al. 1 12009). indirectly validating such a component 
decomposition. 

With the measured velocity Vb.75 = 51.3±4.1 km s _1 of 
C iv absorption lines tracing the BH motion and considering 
the systemic velocity (Vs = +310 km s" 1 ) of LMC X-3 
and velocity semi-amplitude of the companion star Kc = 
235 km s~ , we can constrain the companion to black hole 
mass ratio. At phase 0.75, the BH is receding from us at the 
maximum speed, and thus the relation of V0.75 and Kbh can 
be expressed as 

Vo.75 = Vwmd + Vs + Kbh , (4) 

where V w i n d is the intrisinc velocity of the accretion disk 
wind in our LOS. Assuming Vwmd = —400 km s~\ a similar 
wind velocity found in GRO J 1655-40 (Miller et al. 2006), 
from Eq. [4] we obtain Kbh ~ 140 km s _1 . Plugging these 
numbers in Eq. [2] we further obtain q = Mc/Mbh ~ 0.6, 
which is consistent with the ratio usually adopted for the 
LMC X-3 system (e.g., C83). We therefore suggest that the 
observed C iv absorption feature is consistent to our model 
that the accretion disk wind moves with the BH and the 
observed Doppler shift is a combination of the wind velocity 
and the BH's orbital velocity in our LOS. 
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In the above exercise, we assumed a V w i n d in order to 
obtain Kbh, since there is a degeneracy of Vwind and -Kbh 
in Eq.H Unfortunately, the HST/COS G160M observation 
analyzed above only lasted for one HST orbit, and thus 
cannot break the degeneracy, i.e., we cannot put an inde- 
pendent measure of the BH's velocity, without observing 
the orbital modulation of the Doppler motion of the wind. 
The HST/COS G130M observations in a shorter wavelength 
band covered a much larger part of the orbital period. As 
shown in Fig. [4] the combined spectrum revealed several 
complicated absorption features. Nevertheless, Si n and Si iv 
absorption lines at Vb.75 ~ 50 km s _1 are also detected; 
the apparent higher velocity components may arise from 
the outer part of the accretion disk with lower local veloci- 
ties, mimicking the different velocities a nd ionization zone s 
of AGN warm-absorbers/outflows (e.g., lArav et al.ll2005T ). 
However, the combination of the complexity of the observed 
absorption features and the rather incomplete orbital phase 
coverage of these previous HST/COS observations does not 
warrant further more quantitative analysis for breaking the 
above mentioned degeneracy and probing the nature of those 
higher velocity components. 

Clearly, observations covering more orbital phases are 
badly needed to break the degeneracy. Observations around 
phase 0.25 are the most favorable ones for this purpose. In 
contrast to the existing observations taken around 0.75 in 
which the BH is receding at the maximum velocity from us, 
at phase 0.25, the BH is expected to be moving toward us 
at the maximum velocity, so is the disk wind. Therefore, the 
relation between V0.25 and Kbh can be expressed as 

V0.25 = V„md + V S - Kbh- (5) 

If Vs = —400 km s _1 as assumed, the Vb.25 is expected to 
be at —230 km s _1 . The real measurement of Vb.25 from the 
future observations, together with Vb.75 measured from the 
existing observations, would allow us to solve Vwind and Kbh 
from Eqs. [4] and [5] and then to reliably constrain Mc/Mbh 
and system inclination angle i (Eqs. [2]and[T]). 



5 SUMMARY AND DISCUSSION 

As shown in Eq. [2 Mbh can be directly obtained if Kbh 
can be measured, in addition to Mc and Kc measured by 
observing the companion star. In this work we suggest to 
measure Kbh by detecting the Doppler orbital motion of 
the accretion disk wind absorption lines, assuming that the 
accretion disk wind moves with the BH and does not have 
systematic orbital phase dependence. This method has the 
potential of circumventing the model dependence and other 
uncertainties in estimating the orbital inclination angle, that 
is required in the method commonly used to measure the BH 
masses in XRBs by detecting the companion star's velocity 
and light curves, as shown in Eq. [1] Actually knowing the 
mass ratio with Eq.[5] one can in turn use Eq.[T]to derive the 
inclination angle, which can be used to calibrate the light 
curve model used previously to derive the inclination angle. 

Our analysis of the previous Chandra/HETG observa- 
tions of GRO J 1655-40 have revealed wind velocity mod- 
ulation consistent with the orbital motion of the BH pre- 
dicted from its previously measured system parameters. An 



independent projected radial velocity semi-amplitude mea- 
sured here allows its inclination angle determined without 
using the modeling of its ellipsoidal light modulation of its 
companion. We find its BH radial velocity semi-amplitude 
Kbh = 90.8 ± 11.3 km s" 1 , BH mass Mbh = 5.41±g;gf M© 
and system inclination i — 72.0 ° , where Mbh does not 
depend on i at all. However the very limited orbital coverage 
of the observations does not allow more accurate system pa- 
rameter measurements of this binary system. Nevertheless 
with the velocity component of its orbital motion removed, 
we can obtain more accurate measurements of the intrinsic 
velocities of each line along our line of sight, and thus may 
be able to constrain further the physical properties of the 
wind, by combining with the velocity broadenings of these 
lines; this is the subject our future work. 

Our analysis of the previous HST/COS observations of 
the HMXB LMC X-3 has found absorption line features 
consistent with that predicted by assuming the previously 
measured dynamical parameters of LMC X-3 and the wind 
properties in LMC X-3 being similar to that observed in 
another BHXB GRO J1655-40. Given the limitations of the 
previous HST/COS observations of LMC X-3 that do not al- 
low to break the degeneracy of the wind velocity and the BH 
orbital velocity, new HST/COS observations are required to 
cover significantly different orbital phases. 

As mentioned in S10, the C iv features shown in Fig. [3] 
might be P-cygni profiles. However, the emission features 
agree well with all other emission lines detected, which are 
attributed to the heated stellar surface. Therefore it is more 
reasonable to attribute the emission features in Fig.[3]to the 
heated stellar surface (as done in S10), thus invalidating the 
P-cygni profile interpretation. UV emission lines produced 
from the heated stellar surface are not unique in the sys- 
tem of LMC X-3, but rather a common feature observed in 
XRBs (e.g., Vrtilek et al. 2003). Of course new observations 
suggested above would definitely reveal the nature of the the 
C iv features shown in Fig. [3] 

In this work, we have also assumed that the accretion 
disk wind velocity is constant, at least during one full or- 
bital phase, in order to apply this method reliably. In re- 
ality, the intrinsic wind velocity may have random fluctu- 
ations, though the fluct uations do not se em to be signifi- 
cant in GRO J1655-40 |Miller et alj|2006h . However, it has 
been known that wind absorption features are not always 
detected and it is also not fully understood when and why 
wind absorptions are present or absent. Future high signal 
to noise observations covering more orbital phases may shed 
some lights on this problem and test ultimately if our sug- 
gested method can be applied reliably and produce accurate 
BH mass measurements in BHXBs. The joint JAXA/NASA 
ASTRO-H mission is particularly suitable for making such 
observations, with its high-throughput spectroscopy pro- 
vided by the mic ro-calorimeter with hig h spectral resolution 
of AE ~ 7 eV dTakahashi et alj|201oh . Finally we should 
point out that this method in principle can also be applied 
to other accreting compact objects, such accreting neutron 
star and white dwarf binaries, with detectable accretion disk 
wind absorption line features. 
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ABSTRACT 

So far essentially all black hole masses in X-ray binaries have been obtained by 
observing the companion star's velocity and light curves as functions of the orbital 
phase. However a major uncertainty is the estimate of the orbital inclination angle of 
an X-ray binary. Here we suggest to measure the black hole mass in an X-ray binary 
by measuring directly the black hole's orbital motion, thus obtaining the companion 
to black hole mass ratio. In this method we assume that accretion disk wind moves 
with the black hole and thus the black hole's orbital motion can be obtained from 
the Doppler velocity of the absorption lines produced in the accretion disk wind. We 
validate this method by analyzing the Chandra/HETG observations of GRO J1655-40, 
in which the black hole orbital motion (Kbh = 90.8 ± 11.3 km s _1 ) inferred from the 
Doppler velocity of disk-wind absorption lines is consistent with the prediction from 
its previously measured system parameters. We thus estimate its black hole mass 
(M BH = 5.41±oj7 M©) and then its system inclination (i = 72.0t 7 7 8 5 °), where M BH 
does not depend on i. Additional observations of this source covering more orbital 
phases can improve estimates on its system parameters substantially. We then apply 
the method to the black hole X-ray binary LMC X-3 observed with HST/COS near 
orbital phase 0.75. We find that the disk-wind absorption lines of C iv doublet were 
shifted to ~ 50 km s _1 , which yields a companion-to-black-hole mass ratio of 0.6 
for an assumed disk wind velocity of —400 km s _1 . Additional observations covering 
other orbital phases (0.25 in particular) are crucial to ease this assumption and then 
to directly constrain the mass ratio. This method in principle can also be applied 
to any accreting compact objects with detectable accretion disk wind absorption line 
features. 



1 INTRODUCTION 

Black holes (BHs) are believed to exist in many X-ray bi- 
naries (XRBs) and active galactic nuclei (AGN). Measure- 
ment of the motion of a BH with respect to its surrounding 
in such a system can in principle set strong constraints on 
the mass of the BH. However in an AGN this is normally 
not possible because the time scale of BH's significant mo- 
tion is too long and/or the BH is barely moving at all. In 
an XRB, the BH should move with respect to the system's 
center of mass (CM), making it possible to detect directly 
the BH's motion with respect to the CM. However in a low- 
mass XRB (LMXB) the companion is normally much less 
massive than the BH, so that the BH barely moves or moves 
very slowly with respect to the CM. In some LMXBs, the 
companion's mass is comparable or only sever al times less 
mass i ve than the BH, e.g., in GRO J1655-40 (|Zhang et al.l 
1 19941 ; lOrosz fc Bailvnlll997lh the BHs' motion may be sig- 



nificant enough for direct detection. The most favorable sys- 
tems for detecting BH's motion should be high-mass XRBs 
(HMXBs), in which the BHs move rapidly with respect to 
their CMs. 

Orbital motion of double-peaked disk emission 
line s were observed for n eutron star XRB Sco X- 
1 ilSteeghs fc Casaxesl l2002r), the BH X RB A0620-00 
l|Haswell fc Shafterl Il990l: lOrosz et all Il994l ). and the BH 
GRS 1124-68 |Orosz et al.lll994h .~Unfortunatelv a signif- 
icant phase offset of velocity modulation was found from 
that expected based on the observed orbital motion of the 
companion, though the veloci ty semi-amplitude is consistent 
with the expected mass ratio (|Orosz et al.lll994h . Soria et al. 
(1998) observed the orbital motion of the double-peaked disk 
emission line He n A4686 from GRO J1655-40, and found 
its velocity modulation phase and semi-amplitude in agree- 
ment with the kinematic and dynamical parameters of the 
system. Therefore a more robust mass lower limit is placed 
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based on the observed motion of the primary and thus ruling 
out any possibility for a neutron star as the primary in the 
system (Soria et al. 1998). However one major problem in 
accurately measuring the orbital motion of the primary from 
the observed double-peaked emission lines is how to deter- 
mine reliably the line center, because the lines are typically 
asymmetric and also variable. 

In an XRB, both the accretion disk and its wind move 
with the BH, and thus provides us with another opportunity 
to measure the BH's motion via Doppler shift of absorption 
features of the accretion disk wind. Accretion disk winds are 
ubiquitous in XRBs and normally detected through ionized 
absor ption lines, typically with around ~ 10 3 km s -1 o r less 
(e.g.. lUeda etld1l2004l ; [Miller et alJliooj , 120061 . liooj ). but 



can reach to about 



10 km s in some extreme cases, 



e.g. i n the newly dis covered BH transient IGR J17091- 
3624 l|King et alj 201 ll ) . In particular the high quality Chan- 
dra/HETG observations of the BHXB GRO J 1655-40 have 
found many highly ionized narrow absorption lines, which 
are interpreted as evidence of magnetic field-driven ac- 
cretion disk wind dMiller et all 120061 ; iKallman et"aH 120091 ; 
iLuketic et al.|[2010l ); although the absorption lines are also 
interpreted as from the absorption by X-ray heated thermal 
wind (Netzer 2006). Regardless the origin of the accretion 
disk wind, the orbital motion of its many absorption lines 
may be measured reliably, because typically many narrow 
absorption lines are present with high signal to noise ratios 
and appear to be rather stable when observed. 

In this work, we model the Doppler motion of wind ab- 
sorption lines from the LMXB GRO J1655-40 and HMXB 
LMC X-3 to constrain directly the companion to primary 
mass ratio, in order to measure their BH masses and or- 
bital inclination angles. We first describe the methodology 
and test its feasibility by applying it to Chandra/HETG 
observations of GRO J1655-40, revealing for the first time 
the velocity modulation of wind absorption lines in an XRB 
and thus providing a new measurement of its BH motion 
and orbital inclination angle. We then apply this method 
to the HST/COS observations of LMC X-3 attempting to 
constrain the companion to primary mass ratio. Finally we 
discuss further observations needed to achieve the required 
accuracy of BH mass estimate for LMC X-3, as well as po- 
tential problems and uncertainties of applying this method. 



2 METHODOLOGY 

So far, all BH masses in XRBs have been estimated using 
the Kepler's 3rd law of stellar motion, expressed in the so 
called the mass function, 



f(M) EE P olh K c /27vG 



M BH sin 3 i/(l + g) 2 , (1) 



where P or b is the orbital period, Kc is the semi-amplitude of 
the velocity curve of the companion star, Mbh is BH mass, 
i is the the orbital inclination angle, and q = Mc/Mbh 
is the mass ratio. Since the only direct observables are P OT b 
and Kc, both Mc and i have to be determined indirectly in 
order to obtain the BH mass estimate reliably. The compan- 
ion's mass Mc can be determined relatively reliably by the 
observed spectral type of the companion star. For LMXBs, 
i can be estimated by modeling the optical or infrared light 



curve modulation, though model dependence and other un- 
certainties (such as accretion disk contamination) cannot be 
circumvented completely. For HXMBs, i is normally not de- 
termined very well; in many cases observations or lack of 
eclipse of the accretion disk emission by the companion is 
used to put some constraints on the possible ranges of i. For 
de tails of BH mass estimate s using this method, please refer 
to lRemillard fe McClintockl i|2006l i. 

On the other hand, the mass ratio q can be determined 
directly according to the law of momentum conservation, 



Mc/Mbh = Kbb/K c , 



(2) 



if the semi-amplitude of the velocity curve of the BH Kbh 
can be observed directly. Since a BH is not directly observ- 
able, we can only hope to observe any emission or absorption 
line feature co-moving with it. The accretion disk certainly 
moves with the accreting BH. However any line feature of 
the inner accretion disk suffers from the broadening of disk's 
orbital motion and distortions by relativistic effects around 
the BH, thus making it practically impossible, or difficult to 
the least, for detecting the binary orbital motion of the BH. 
Emission line features from the outer disk region are nor- 
mally detected with double-peaks, which can be modeled to 
obtain the semi-amplitude of the velocity of the compact 
object, as discussed above. However a major uncertainty is 
to determine the mean separation between the emission re- 
gions of the blue-shifted and red-shifted components. 

Fortunately, as we have discussed above, absorption line 
features of accretion disk winds in BHXBs have been rou- 
tinely detected with high significance. This suggests that 
the accretion disk wind in an XRB moves with the disk that 
produces the wind, since otherwise the wind from the disk 
would not intercept the continuum emissions produced from 
the same disk. In this case, the Doppler motion of wind 
absorption line features can be considered as that of the 
BH, unless the wind interacts strongly with the surrounding 
interstellar medium or the wind intrinsic velocity has sys- 
tematic orbital dependence. It has been found recently that 
large scale (pc) cavities exist around microquasars (XRBs 
producing relativistic jets), and perh aps are ubiquitous i n 
all BHXBs producing strong winds (|Hao fc Zhand |2009h . 
This suggests that at least in the vicinity of the BH the 
wind cannot interact directly with interstellar medium. The 
interactions of disk winds with interstellar medium at pc 
scale would not complicate the observations, because it takes 
much longer than an orbital period before the winds could 
arrive at the boundaries of these cavities, and thus the wind 
should have lost any possible memories of the orbital motion 
of the compact object. For disk-fed accretion, the wind ve- 
locity is not expected to have any orbital phase modulation; 
for stellar-wind accretion the focused wind may complicate 
the situation considerably though. We therefore suggest to 
apply this method to only those systems in which the accre- 
tion is disk-fed, i.e., no significant wind interaction happens 
between the stellar and accretion disk winds. Nevertheless 
only observations can tell us in which kinds of systems this 
method can be used reliably. 
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Figure 1. Distributions of the line velocities in four equally 
spaced intervals and their average values, during the Chan- 
dra/HETG observations lasting for more than 60 ksec; note that 
here the system velocity of the BH XRB has already been sub- 
tracted. Panels from (a) to (d) corresponding to the phases shown 
in Fig.[2]from left to right. 
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Figure 2. Velocity curve of the 39 observed absorption lines, 
after rejecting six lines as suspected outliers and subtracting the 
line of sight intrinsic velocity of each line at each orbital phase. 
The upper panel marks the velocity of each line with its 1-rr error 
bar slightly shifted horizontally for visual clarity; the inset shows 
all velocities, including the three data points out of the range in 
the main panel. The bottom panel shows the weighted average 
velocity of all lines in the upper panel at each phase; the solid 
curve is the fitted velocity curve with its orbital period and phase 
fixed at the values observed previously. 



3 APPLICATION TO THE LMXB 
GRO J 1655-40: A TEST STUDY 

GRO J1655-40, discovered as an X-ray transient by 
Zhang et al. (1994), is a well-known and best studied 
BHXB and the se cond source with superluminal relativis- 
tic jets detected dTingav et al.l 1 19951 ; lHarmon et al. I Il995l ; 
iHiellming fc Rupenl 1 19951 ). Its system parameters remain 
so far the best measured among all known BHXBs, with 
P orb = 2.62191 ± 0.00020 d, 1/q = M BH /M C = 2.6 ± 0.3, 
70.2±1.9°, and M B h = 6.3±0.5 Mp (al l 95% confidence) 



i|Orosz fc Bailvnlll997l ; iGreene et al.ll200ll ). Its precise BH 
mass and inclination measurements allow its BH spin pa- 
rameter determined from its X-ray continuum fitting, first 
proposed by Zhang et al. (1997) and then refine d by incor- 
porat i ng detailed modeling of various effects (e.g..lYao et al 
20051; iLi et "ail 120051; iDavis et alj 120051; IShafee et al.l I200T 



Steiner et al.l 120091 ; (McClintock et al.l l201l[ ). However pos- 



sible sources of systematic errors are in the modeling of the 
ellipsoidal light modulation of the companion star and con- 
tamination of the disk's continuum emission in the optical 
to infrared bands. Soria et al. (1998) observed the velocity 
curve of the double-peaked disk emission line He 11 A4686 
from GRO J1655-40, and derived its projected radial veloc- 
ity semi-amplitude determined for the primary as Kbh = 
76.2 ± 7.5 km s~\ yielding M BH = 6.62 ± 0.74 M© and 
i = 66.6 ± 7.7° (all 68.3% confidence), fully consistent with 
that determined from the companion's velocity curve and 
ellipsoidal light variation. 

Miller et al. (2006) detected many narrow and ionized 
absorption lines with High Energy Transmission Grating 
(HETG) Chandra observations of GRO J1655-40. The exact 
position of the point source cannot be determined accurately 
in the zeroth order image, which is severely piled-up due to 
its high flux. An offset in the source position will cause a 
systematic shift of all lines in the wavelength space. How- 
ever, this shift is mostly canceled out in the combined total 
spectrum if each line is detected in the spectra of both the 
+ 1 and —1 orders with similar count spectra. However the 
combined lines will be broadened by the unknown system- 
atic offset, which in turn will cause larger uncertainties in 
determining the center of each line. There may also exist 
other sources of systematic errors which can cause offset be- 
tween the centers of the same line in the spectra of the +1 
and —1 orders. 

Here we determine the position of the zeroth order im- 
age by minimizing the relative offsets between the absorp- 
tion lines in the +1 and —1 orders, and estimate the sys- 
tematic error in determining the center wavelength of each 
absorption line. We first fit each absorption line feature with 
a Gaussian profile in the +1 and —1 orders of the High En- 
ergy Grating (HEG) and Medium Energy Grating (MEG) 
separately, i.e., we get a pair of profiles for each absorption 
line, integrated over the whole observation in order to get 
the highest signal to noise ratio. The local continuum around 
each line is assumed to have a power-law shape, but our re- 
sults are insensitive to the continuum shape. For HEG and 
MEG we obtain 43 and 9 pairs of absorption lines with at 
least 3-cr detection, resulting in 43 and 9 differences of the 
central wavelengths of these pairs, respectively. We then use 
the Bayesian method to find the offset d and its systematic 



4 S. N. Zhang, J. Liao & Y. Yao 



error ct s by maximizing the following function, 



p(d,a B ) = 



1 



n 



(3) 



where Xk and o Xk are the wavelength difference and its sta- 
tistical error of the fc-th pair (k = 1 to n, and n = 43 for 
HEG or 9 for MEG), respectively. Finally for HEG we get 
d = (4.8 ± 0.1) x 1CT 3 A and a s = (3.2 ± 1.3) x 1CT 4 A. 

We thus take the error for d as ad — cr B — 3.2 x 1CP 4 A 
for HEG. Similarly for MEG we get d = (9.7 ± 0.8) x 10~ 3 
A for MEG, where <7<j = 9 x 10~ 4 A is dominated by the 
systematic error. Taking these additional errors for MEG 
and HEG line wavelengths and requiring all pairs in HEG 
with d = 4.8 x 10" 3 A and in MEG d = 9.7 x 10" 3 A, 
we have a total \ 2 ~ 55 for 50 degrees of freedom, thus 
validating our method of determining the systematic errors. 
Since the wavelength bin (an image pixel) is 2.5 x 10~ 3 A or 
5.0 x 10~ 3 A for HEG or MEG, we thus shift the spectra of 
the —1 and the +1 orders by one bin, towards the decreasing 
(— 1 order) and increasing (+1 order) wavelength directions, 
respectively. This is identical to shifting the location of the 
point source by one pixel in the zeroth order image. Finally 
we combine the two shifted spectra as one single spectrum. 

Because the observations lasted for about a quarter of 
the orbital phase of GRO J1655-40 and many absorption 
lines have high signal to noise ratios, here we divide the ob- 
servations into four equal orbital internals in order to detect 
the orbital motion of the absorption lines. For each com- 
bined spectrum we again fit each absorption feature with 
a Gaussian profile to determine its central wavelength, line 
width and intensity. Table 1 lists the 45 lines in all four 
intervals with at least 3-<x detection that are included for 
further analysis; here only statistical errors are shown. 

Fig.Q]shows the velocity distributions of the four groups 
of 45 lines and their average velocities; note that here the 
system velocity of the BH XRB has already been sub- 
tracted. Adopting the orbital period and phase ephemeras 
from Greene et al. (2001) and assuming that the BH's mo- 
tion is exactly anti-phased with the companion's motion 
(Soria et al. 1998), we fit velocities of these lines simul- 
taneously to a sinusoidal function. In the fit, the intrinsic 
velocity of each line is a free parameter, but all lines in 
the same orbital phase follow the same orbital modulation. 
This way we have 46 free parameters (45 intrinsic velocities 
for all these lines plus the line of sight velocity amplitude 
of the black hole) with 180 data points. The fit results in 
Kbh = 93.8 ± 11.1 km s _1 with x 2 = 179 for 134 degrees 
of freedom; the systematic errors determined above are in- 
cluded for all lines. 

Considering that some lines may be outliers, we reject 
the four groups of six lines with central velocities more than 
300 km s _1 from the median values of each distribution, 
as listed in Table 1. The remaining four groups of 39 lines 
have velocity dispersion (1-cr) of about 60 km s _1 , i.e., the 
rejected lines are more than 5-<r away from the median val- 
ues. The fit to the remaining four groups of 39 lines yields 
Kbh = 90.8 ± 11.3 km s _1 with \ 2 = 131 for 116 degrees of 
freedom, a marginal improvement over the full dataset fit. 
This means that the deviations of these possible outliers are 
not significant statistically and thus do not deserve further 
studies at this stage. Fig. [2] shows the fitting results; please 
note that the velocities are obtained by subtracting the fitted 
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Figure 3. Spectrum of the available HST/COS observations 
around C IV doublet. Red curves indicate the fit to the emis- 
sion from the heated surface of the companion star, green curves 
are the fit to the Galactic ISM absorptions at km s —1 , and blue 
curves indicate the fit to the accretion disk wind absorptions at 
50 km s~ x (local velocity S3 —400 km s _1 ). The vertical bars 
mark the central positions of these components. 



line of sight intrinsic velocity of each line from the obtained 
central velocity of each line at this orbital phase. It is worth 
noting that the phase zero in this system was defined as that 
when the companion is receding from the observer at the 
maximum velocity (Orosz & Bailyn 1997). Nevertheless both 
results are statistically consistent with Kbh = 76.2 ±7.5 km 
s" 1 obtained by Soria et al. (1998). If we fit the four groups 
of 39 lines to a straight line, we get \ 2 = 813 for 154 de- 
grees of freedom. The linear model is thus rejected with high 
significance, compared to the sinusoidal model. 

With the fitted Kbh = 90.8 ± 11.3 km s" 1 and tak- 
ing the other system parameters (except the inclination) of 
GRO J1655-40 from Greene et al. (2001), we first obtain 
its BH mass M BH = 5Al±°- 9 £ M Q from Eq. 2 and then 
its system inclination i — 72.OI75 from Eq. [1] respec- 
tively. All these (l-o") errors are obtained by Monte-Carlo 
samplings, because of the asymmetry and coupling of some 
errors. These parameters, albeit with large uncertainties due 
to the very incomplete orbital coverage of observations, are 
consistent with all previous measurements. Therefore the 
existing data of GRO J 1655-40 validate our proposal that 
absorption lines produced in the accretion disk wind can be 
used to measure directly the orbital motion of the BHs in 
BHXBs. 



4 APPLICATION TO THE HMXB LMC X-3 

LMC X-3 is another excellent object for applying this 
method. It is a bright XRB system in the Large Magellanic 
Cloud composed of a B3 V star and a central BH, and reced- 
ing away from us at a systemic velocity Vs = +310 km s _1 
l|Cowlev et al. I Il983l ; hereafter C83). It is one of the few 
BH systems that are persistently luminous in both X-ray 
and far-ultraviolet (FUV) wavelength bands. The BH is be- 
lieved to be undergoing accretion from its B-star compan- 
ion via Roche lobe overflow with an orbital period of 1.7 
days (C83); we therefore do not expect any significant con- 
tamination of the stellar wind to any absorption features 
of its accretion disk wind. Spectroscopic observations of 
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Figure 4. Si II and Si IV absorption features in the spectrum 
of the previous HST/COS G130M observations, normalized to 
the continuum flux. The apparent higher velocity components 
(marked with pink curves) may be produced in the outer part of 
the disk wind at lower local velocities. 

the B star indicate a large radial velocity semi-amplitude, 
Kc = 235 ± 11 km s" 1 (C83). Taken the mass of the B3 V 
star as about 6 Mq, the BH mass in LMC X-3 is thus es- 
timated to be 5 — 10 Mp,, assuming an inclination an gle of 
50° - 70° (C83; iKuiper et all [l98s[ ISoria et al]|200ll ). Be- 
cause of the considerable uncertainty in its inclination (and 
thus BH mass), its BH spin has not been reliably deter- 
mined yet with the X-ray continuum m odeling method (e.g., 
IZhang et aflll997l : lsteiner et alJl20ld ). Since the BH should 
be moving at a comparable velocity to its companion, it is 
possible to directly measure its orbital motion if the absorp- 
tion lines produced in its accretion disk wind are detected. 

Although LMC X-3 has been intensively studied in the 
X-ray band, the previously operated and currently operat- 
ing X-ray spectroscopic instruments lack the combination 
of the required sensitivity and resolution to measure the 
expected Doppler motion of any accretion wind absorption 
features. High resolution UV spectroscopic observations of 
this source are rather sparse; so far only two FUSE ob- 
servations and four observations with the Cosmic Origins 
Spectrograph (COS) onboard the Hubble Space Telescope 
(HST) (HST/COS her e after) to LMC X-3 have been made 
jHutchings et al.ll2003l ; IWang et al.ll2005l ; ISong et allkoid . 
here after S10). Again FUSE does not have the combination 
of the required sensitivity and the observations did not pro- 
duce spectra with high enough signal to noise ratio for the 
purpose of this investigation. The HST/COS does have the 
required performance to do so, although the available obser- 
vations only covered a small portion of the orbital period of 
LMC X-3. In this work, we analyze the available HST/COS 
UV spectroscopic observations of LMC X-3, aiming at con- 
straining the systemic parameters and further demonstrat- 
ing the feasibility of our proposed new method of measuring 
the BH's mass in an XRB. 

From the FUSE and HST/COS observations, S10 de- 
tected variable O VI and N v emission lines and found that 
variability of their intensities are inconsistent with expecta- 
tion of a stellar wind origin. So they attributed them to the 
heated stellar surface of the companion star. In fact both 
absorption and emission features are detected in N v (with 
G130M) and C iv (with G160M) doublets (Fig. 2 in S10) in 
the HST/COS observations (Program 11642; please refer to 



S10 for detailed description of the HST/COS observations 
and data analysis). However, since S10 focused on the ori- 
gin of the emission variability and the connection between 
the emission features and the systemic parameters, these 
absorption features have not been explored in detail. 

Here we analyze the C iv doublet complexes (we do not 
conduct the same analysis to the N v features because the 
absorptions are much less significant than those of C iv). 
In Fig. [31 we decompose the observed C iv doublet into 
three components for each line: one emission component and 
two absorption components. The emission component is as- 
sumed to come from the heated surface of the companion, 
the same as the O VI and N v emission lines (S10); its width 
and Doppler shift are fixed to those inferred from the O VI 
and N v emission lines in the decomposition. One absorp- 
tion component (at zero velocity) is assumed to come from 
the Galactic ISM absorption. The other significant absorp- 
tion component has a redshift of about 50 km s _1 , and may 
have two possible origins. One is the absorption produced in 
the stellar wind of the B3 V companion star. However, the 
orbital phase of the binary system was about 0.75 during 
the HST/COS observations, i.e., the straight line connect- 
ing the companion star and the BH is perpendicular to our 
LOS (please be noted that, in contrast to GRO J1655-40, 
the orbital phase zero in LMC X-3 was defined as that when 
the BH is at its superior conjunction, i.e., the companion is 
just between the observer and the BH; Cowley 1983; S10). 
In order for any stellar wind to intercept significantly the 
emission from the inner disk region, a significant amount 
of stellar wind has to stream to the inner disk region, i.e., 
the system must be wind-fed, against the common believe 
that LMC X-3 is actually disk-fed. We thus rule out this 
possibility. The other and the only viable scenario is the ab- 
sorption by the accretion disk wind, making it possible to 
directly measure the orbital motion of the BH. It is worth 
noting that the measured C iv column density in the ISM 
[log A r isM(cm -2 ) = 13.7 ± 0.1] agrees remarkably well with 
predicted value [log JVciv(cm -2 ) = 13.8] in a joint analysis 
of the X-ray and FU V spectroscopic observations of LMC X- 
3 (|Yao et al. 1 12009). indirectly validating such a component 
decomposition. 

With the measured velocity Vb.75 = 51.3±4.1 km s _1 of 
C iv absorption lines tracing the BH motion and considering 
the systemic velocity (Vs = +310 km s" 1 ) of LMC X-3 
and velocity semi-amplitude of the companion star Kc = 
235 km s~ , we can constrain the companion to black hole 
mass ratio. At phase 0.75, the BH is receding from us at the 
maximum speed, and thus the relation of V0.75 and Kbh can 
be expressed as 

Vo.75 = Vwmd + Vs + Kbh , (4) 

where V w i n d is the intrisinc velocity of the accretion disk 
wind in our LOS. Assuming Vwmd = —400 km s~\ a similar 
wind velocity found in GRO J 1655-40 (Miller et al. 2006), 
from Eq. [4] we obtain Kbh ~ 140 km s _1 . Plugging these 
numbers in Eq. [2] we further obtain q = Mc/Mbh ~ 0.6, 
which is consistent with the ratio usually adopted for the 
LMC X-3 system (e.g., C83). We therefore suggest that the 
observed C iv absorption feature is consistent to our model 
that the accretion disk wind moves with the BH and the 
observed Doppler shift is a combination of the wind velocity 
and the BH's orbital velocity in our LOS. 
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In the above exercise, we assumed a V w i n d in order to 
obtain Kbh, since there is a degeneracy of Vwind and -Kbh 
in Eq.H Unfortunately, the HST/COS G160M observation 
analyzed above only lasted for one HST orbit, and thus 
cannot break the degeneracy, i.e., we cannot put an inde- 
pendent measure of the BH's velocity, without observing 
the orbital modulation of the Doppler motion of the wind. 
The HST/COS G130M observations in a shorter wavelength 
band covered a much larger part of the orbital period. As 
shown in Fig. [4] the combined spectrum revealed several 
complicated absorption features. Nevertheless, Si n and Si iv 
absorption lines at Vb.75 ~ 50 km s _1 are also detected; 
the apparent higher velocity components may arise from 
the outer part of the accretion disk with lower local veloci- 
ties, mimicking the different velocities a nd ionization zone s 
of AGN warm-absorbers/outflows (e.g., lArav et al.ll2005T ). 
However, the combination of the complexity of the observed 
absorption features and the rather incomplete orbital phase 
coverage of these previous HST/COS observations does not 
warrant further more quantitative analysis for breaking the 
above mentioned degeneracy and probing the nature of those 
higher velocity components. 

Clearly, observations covering more orbital phases are 
badly needed to break the degeneracy. Observations around 
phase 0.25 are the most favorable ones for this purpose. In 
contrast to the existing observations taken around 0.75 in 
which the BH is receding at the maximum velocity from us, 
at phase 0.25, the BH is expected to be moving toward us 
at the maximum velocity, so is the disk wind. Therefore, the 
relation between V0.25 and Kbh can be expressed as 

V0.25 = V„md + V S - Kbh- (5) 

If Vs = —400 km s _1 as assumed, the Vb.25 is expected to 
be at —230 km s _1 . The real measurement of Vb.25 from the 
future observations, together with Vb.75 measured from the 
existing observations, would allow us to solve Vwind and Kbh 
from Eqs. [4] and [5] and then to reliably constrain Mc/Mbh 
and system inclination angle i (Eqs. [2]and[T]). 



5 SUMMARY AND DISCUSSION 

As shown in Eq. [2 Mbh can be directly obtained if Kbh 
can be measured, in addition to Mc and Kc measured by 
observing the companion star. In this work we suggest to 
measure Kbh by detecting the Doppler orbital motion of 
the accretion disk wind absorption lines, assuming that the 
accretion disk wind moves with the BH and does not have 
systematic orbital phase dependence. This method has the 
potential of circumventing the model dependence and other 
uncertainties in estimating the orbital inclination angle, that 
is required in the method commonly used to measure the BH 
masses in XRBs by detecting the companion star's velocity 
and light curves, as shown in Eq. [1] Actually knowing the 
mass ratio with Eq.[5] one can in turn use Eq.[T]to derive the 
inclination angle, which can be used to calibrate the light 
curve model used previously to derive the inclination angle. 

Our analysis of the previous Chandra/HETG observa- 
tions of GRO J 1655-40 have revealed wind velocity mod- 
ulation consistent with the orbital motion of the BH pre- 
dicted from its previously measured system parameters. An 



independent projected radial velocity semi-amplitude mea- 
sured here allows its inclination angle determined without 
using the modeling of its ellipsoidal light modulation of its 
companion. We find its BH radial velocity semi-amplitude 
Kbh = 90.8 ± 11.3 km s" 1 , BH mass Mbh = 5.41±g;gf M© 
and system inclination i — 72.0 ° , where Mbh does not 
depend on i at all. However the very limited orbital coverage 
of the observations does not allow more accurate system pa- 
rameter measurements of this binary system. Nevertheless 
with the velocity component of its orbital motion removed, 
we can obtain more accurate measurements of the intrinsic 
velocities of each line along our line of sight, and thus may 
be able to constrain further the physical properties of the 
wind, by combining with the velocity broadenings of these 
lines; this is the subject our future work. 

Our analysis of the previous HST/COS observations of 
the HMXB LMC X-3 has found absorption line features 
consistent with that predicted by assuming the previously 
measured dynamical parameters of LMC X-3 and the wind 
properties in LMC X-3 being similar to that observed in 
another BHXB GRO J1655-40. Given the limitations of the 
previous HST/COS observations of LMC X-3 that do not al- 
low to break the degeneracy of the wind velocity and the BH 
orbital velocity, new HST/COS observations are required to 
cover significantly different orbital phases. 

As mentioned in S10, the C iv features shown in Fig. [3] 
might be P-cygni profiles. However, the emission features 
agree well with all other emission lines detected, which are 
attributed to the heated stellar surface. Therefore it is more 
reasonable to attribute the emission features in Fig.[3]to the 
heated stellar surface (as done in S10), thus invalidating the 
P-cygni profile interpretation. UV emission lines produced 
from the heated stellar surface are not unique in the sys- 
tem of LMC X-3, but rather a common feature observed in 
XRBs (e.g., Vrtilek et al. 2003). Of course new observations 
suggested above would definitely reveal the nature of the the 
C iv features shown in Fig. [3] 

In this work, we have also assumed that the accretion 
disk wind velocity is constant, at least during one full or- 
bital phase, in order to apply this method reliably. In re- 
ality, the intrinsic wind velocity may have random fluctu- 
ations, though the fluct uations do not se em to be signifi- 
cant in GRO J1655-40 |Miller et alj|2006h . However, it has 
been known that wind absorption features are not always 
detected and it is also not fully understood when and why 
wind absorptions are present or absent. Future high signal 
to noise observations covering more orbital phases may shed 
some lights on this problem and test ultimately if our sug- 
gested method can be applied reliably and produce accurate 
BH mass measurements in BHXBs. The joint JAXA/NASA 
ASTRO-H mission is particularly suitable for making such 
observations, with its high-throughput spectroscopy pro- 
vided by the mic ro-calorimeter with hig h spectral resolution 
of AE ~ 7 eV dTakahashi et alj|201oh . Finally we should 
point out that this method in principle can also be applied 
to other accreting compact objects, such accreting neutron 
star and white dwarf binaries, with detectable accretion disk 
wind absorption line features. 
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A1 1 /9S^ 1 7fi/zL1 ^ 
— 41 1 l^Zo ) 1 1 01^41 ) 


v 


Mg XII Is - 5p 


6.5801 


-510(17) 164(28) 


-502(16) 146(24) 


-510(17) 178(24) 


^91 /90^ 1 01 (1 1 ^ 

— ozi i^zu ) lyi^oij 


v 

1 


Mg XII Is - 3p 


7.1062 


-365(16) 91(31) 


-406(17) 114(25) 


-398(20) 114(32) 


-381(18) 61(49) 


Y 


Mg XII Is - 2p 


8.4210 


-410(9) 339(11) 


-427(9) 278(10) 


-467(9) 259(10) 


-471(10) 263(12) 


Y 


Ne X Is - 4p 


9.7082 


-359(25) 151(46) 


-390(21) 107(39) 


-439(19) 95(30) 


-399(21) 141(32) 


Y 


Na XI Is - 2p 


10.0250 


-346(19) 166(26) 


-339(19) 146(25) 


-443(22) 151(27) 


-441(20) 127(26) 


Y 


Nc X ls7p 


9.2912 


-418(63) 84(93) 


-447(52) 93(84) 


-380(74) 199(63) 


-391(46) 100(80) 


Y 


Nc X Is - 6p 


9.3616 


-776(25) 287(31) 


-846(25) 291(32) 


-837(34) 250(45) 


-945(32) 114(47) 


N 


Ne X Is - 5p 


9.4807 


-365(28) 247(32) 


-322(28) 194(34) 


-370(27) 222(34) 


-391(26) 201(33) 


Y 


Ne X Is - 3p 


10.2389 


-363(18) 222(25) 


-439(20) 229(22) 


-407(15) 143(21) 


-449(16) 207(20) 


Y 


Nc X Is - 2p 


12.1330 


-441(21) 281(25) 


-420(21) 225(25) 


-452(20) 236(26) 


-482(17) 213(27) 


Y 


Ne II 


14.6310 


-389(34) 182(42) 


-393(33) 132(46) 


-410(30) 82(73) 


-457(46) 223(49) 


Y 


O VIII Is - 4p 


15.1762 


-296(90) 361(113) 


-329(81) 314(133) 


-471(68) 186(101) 


-374(42) 41(73) 


Y 


O VIII Is - 3p 


15.9870 


-23(63) 151(90) 


-9(77) 234(79) 


-110(74) 182(111) 


-2(80) 200(94) 


N 



